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Abstract 


We show that for the two-dimensional multiplicative Brun’s algorithm, 
the exponent of convergence is 1 + d, i.e. there is a d > 0 such that for 


almost all x = (xi,^), 


x; -■ 


p) 




rW 


< 


1 


ond Lyapunov exponent is negative. 


(P {,) ) 


1 +d 


(i = 1,2). Thus the sec- 


* 

In 1993, J. C. Lagarias has shown how to use multiplicative ergodic 
theorems to determine the approximation exponent 1 + d for multi¬ 
dimensional continued fractions. Ito, Keane & Ohtsuki 1993 proved 
that for the two-dimensional modified Jacobi-Perron algorithm 
d > 0 (see also Fujita, Ito, Keane & Ohtsuki 1996). In Schweiger 1996, 
a classical result of Paley & Ursell 1930 was used to determine the expo¬ 
nent of convergence of the Jacobi-Perron algorithm in two dimensions. 
Meester 1997 gave another proof for the result on Podsypanin’s modifica¬ 
tion. 

In this paper, we will show that a similar method can be applied to 
Brun’s algorithm in two dimensions. Clearly, this is no surprise since 
Brun’s multiplicative algorithm is a factor of the modified algorithm. 
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We will start with a description of the algorithm; for general references 
on Brun’s algorithm see e.g. Schweiger 1995. 

Definition. Let B = {(xi,x 2 ) 1 > jxti > >r 2 > 0}; Brun’s Algorithm is 
generated by a map S B ( j , N) —> B, where 


S(x u x 2 ) = < 


if— -N>— [j = 1] 

X\ X\J X\ X\ 

(^,— -n) if— -N<^ [j = 2] 

\*i *\ / *1 *1 


N := 


1 


> 1. 


Let := N^'-^xP.x^)) and/(0 :=>(J’ / - 1 (*i 0) >* 2 0) )); 


W 


ify'(/+l) = l,thenx| /+1) = -h - N (/+1) , x^ /+1) = 


(/) ^ 2 (/) 
*1 *1 

W 


VO A 

ify(/+ 1) = 2, thenx| /+1) = ^y,xj> /+1) = iV (/+1) 


x 


x 


M 


The matrices of Brun’s Algorithm are given as follows: 

Definition. Let t > 1; 


then 


and, for t > 2, 


A W - 

IV B * — 


'nm 2 - 7 (/) y(/)-r 


i 




0 


0 


0 >(/)-! 2->(/). 


:=E, 


A (1) 

, (0) 

A ,} 

.(1) 

Ao) 

A (-i) 

At 

At 

At 

U" 

A? 

A2 _1) 


nj } = 


p? p\ n aP 

^ pV pV) 


— Q (/_1) A (/_1) 

\L X 3 


0 ) 
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Hence, for i = 1,2, we get 


; (0 

qW + q( f '} + 


( 2 ) 


Define /* as the largest integer such that /* < /, and./(/*) = 2 (if there is 
no such /* < /, then /* := —1); consequently, (/ + 1)* is defined as the 
largest integer such that (/ + 1)*</ + 1, and j((t + 1) *) =2. Then if 
j[t) = 1, (/ + 1) / in Definition (1) equals/, and (/ + 1) 7/ = /* = (t + 1)*; 
in the other case, we have (/ + 1) / = /*, and (/ + \) n = t = (/ + 1)* 
Hence 

ify(/) = 1 q {t+2) = N (/+1) ? (/+1) + # (/) , and (3) 

if j(t) = 2: ^ (/+2) - N (/+1 V (/+1) + 7 (/ * } (4) 

Of course, (3) and (4) remain valid if we replace q^ ) by p • ^ for i = 1,2. 
Since the following results hold for both p\ ^ and p\ \ from now on 

we will only write p ( ) instead. We continue with a modification of the 
arguments of Paley & Ursell 1930. 


Definition. 


q ( /+1) 

<7 W 

D 1 - 

A +1) 


p lt — 




P [t) 

» 1 /+! 


P (/,) 

> 1 /+! * — 

P w 

P { '*> 


By Eqs. (3) and (4) we get the following relations: 

if At) = 1: P, + 2 = -P, + u P', +2 = N^p' /+ , - P'; +v P" f+2 = P' t+] 

( 5 ) 

if j{t) = 2: P /+2 = -P' f+V P' f+2 = - P'! +v P" t+2 = P /+1 

( 6 ) 

Definition. 

/ip, i ip;i\ m 

f"= m “{jw'TwJ (7) 

Lemma. 

l^/ + i| <? W P/ (8) 


Prcw/: We use (5), (6) and Definition (7): |P/+i| < max{|P,|, |P'|} < 

q {t) Pf 
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Lemma. 

|P^,|<f W P/ (9) 

Proof: Similar to the previous lemma: |P'' +1 1 < max{|P/|, |P'|} 

<? ( V 

Definition. 

B 2 := B fl ((x u x 2 ) :y'(xi,x 2 ) = 2) 

2 

M := 

/'=0 

Let / 0 := min{/ > 0 S f ~'(x\,x 2 ) G M}, /,„+1 tnin{/ > t m 

[x\ , x 2 ) G M}; we thus have y(/ w ) = 2,j(t m + 1) = 2 and 
y(/ w + 2) = 2. Hence, in choosing (x\,x 2 ) G M, we avoid /* being 
too far away from /, which will simplify the following estimates. 

Lemma. 

H(M) > 0 

Proof: We consider the subset M* C M, where for (xi,x 2 ) G M*, 
N(xi,x 2 ) = 1, N(S(x\,X 2 )) = 1 and N(3’ 2 (xi,x 2 )) = 1; clearly, 
/i(M) > jjb(M*). Since all the cylinders J3(y, N) are proper, i.e. 
S(B(j\N)) = B (see e.g. Schweiger 1998), we can apply the local 
inverse 

t / (y-2,N.|)(j'l,J , 2) = 

to the points ( 0 , 0 ), (1, 0 ) and (1.1). We get a triangle whose vertices are given 

and F 3 (l,l) = 


Lemma. 

max {P/„+3, P/,+2, P/,+1} (10) 



/ 1 i \ / *2 

by the points K 3 (0,0) = K 3 (1,0) = (- 

/■x i\ \ z z / V“ 


-, - J , which clearly is of positive measure. 
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Proof: We apply Lemma (8) and Eq. (4) in (1), and similarly relation (6), (8) 
and (4) in (2) 

(1) |P,„ +4 | < 3 < (? (/ '" +4) - 4 {, - ] )p,. + S 

(2) \P' tw+A \ < N<'- +3 >|P, -+3 | + \Pl +3 \ 

< N<'- +3 >|P,„ +3 | + \P fa+2 \ 

< (N (/ '" +3) ^ (/ " +2) +^ (, "' +1) )max{p //ii+ 2 ,p // „ + i} 

< (JV ('-r+3) jV ( / /^+ 2 )q ('-+ 2 ) q- j\j-(/^/+3)^(/ w ) 

max{p //(i+2 ,P/, J ,+i} 

< (n('» + V» +3 > +^" +1 > - q {,m) )maz{pt.+z,Pt.+\} 

< (^ (/ "' +4) -£ ( ''hmax{p,, + 2,p /<(+ i} 

Lemma. 

( 9 (,m) \ 

P /,„+5 < (l ~ (/jt+5) J max{P/.+3, P/,+2, P/,+i} 

Proof: In (1) we use (5), (6) and the previous lemma; (2.1) follows from (5), 
the previous lemma, Lemma (9) and Eq. (3) 

(1) |P,. +5 | < max{|P,,„ +4 |, |P^ +4 |} < (? ( '"' +4) - q [t '" ] ) 

max {P/,.+3, Pt,„+2i Pt,„+ 1} 

(2.1) + 3) = 1: 

|P;„, +5 I < N^ +4) |P^ +4 | + |P" w+4 | 

< (]V (/ - + 'V'- +4) -^('-^maxlp^.p^.p^+i} 
+ <? ('» + V„+3 

< (jv (/ " +4) ^ (/ " +4) + ^ ( '’" +3) - 

max {P/»+3> Pl m +2i Pt,„+ 1} 

< ( ? (/ » +5) - ? (/ » ) )max{p^ + 3,p^ +2 ,p / „ + i} 

(2.2) y(/ w + 3) = 2 Similar to (2) in the previous lemma 
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um.at 


Since Lemma [10] guarantees that pt m +4 < max{p /w + 3 , p /w + 2 , 
for t + 6 we similarly get the following result: 


Lemma. 


P /,+6 < 1 


7 ('-) 


(/,+ 6 ) 




max {P/„+3> Pl„+2, Pt„,+ 1} 


Now let / > t m + 6; we have 

Lemma. 


Pt < 1 - 


7 (/«) 


7 (/,+ 6 ) 


^ max{p /<>(+ 3 , P/.+ 2 , P/„+i} 


Proof : We proceed inductively, where the assumptions are given by the 
previous lemmas; (1) follows from (5), (6) and Definition (7); in (2.1) we 
apply (5), (9) and (3), in (2.2.1) we use (6), (5), Lemma (9), (3) and (4), in 
(2.2.2) (6), Definition (7), Lemma (9) and (4) 


l+l < max{|P / _ 1 |,|P / / _ 1 |} 


< 



A >m] \ 

q{t,„+(>) J 


max {P/„+3) P/.+ 2> P/,+1 } 


(2.1) j(t — 2) = 1: 


|p;i , N<'-'>ip^,i + iPtii 

^ + q ( ' , ~ 2 ^)ma.x{p,-u P/-i} 

qi') 

< maxfp,-!,/)^} 


< 1 




7 (/,+ 6 ) 


max 


{P/»+3. P/„+2, P/,+1 } 


( 2 . 2 ) 


j{t-2) =2: 
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(2.2.1) j{t- 3) = 1 


\K\ N^ip^i + ipy,,! 

a(') - 


< 


< 


< 


< 


qW 

n( / - , )|p / _ 2 | + |p; , _ 1 | 

(N('-Vg('- 3 ) + g('- 2 ))max{p,-2, Pt- 3} 
{N {t ^N (, - 2 ) q (> - 2) + N( / - 1 V (/_3) )max{p / _ 2 ,p / _ 3 } 
jV('-i)^('-i) ma x{p/_2, p t - 3} 


< maxl/J^p,^} 


7 (C//) 


< 1 - 


p('#+6) 


max {P/,„+3> Pt,„+2i Pt,„+ 1 } 


(2.2.2) j(t-3) = 2: 


l-P/l < N( /-1 )|P/_! | + JP^-ij 


/W 


< 


< 


?w 


( /—1 ) |p /—1 1 + |P/_ 2 | 


7 W 


(N (/ 1) ^ (/ 3) )max{ < o / _ 1 ,p / _ 3 } 




W 


< max{p,_i, p,_ 3 } 

] max{p /)(;+3 , pt m +ii Pt„,+\} 


< 1 


r 


7 (/ w +6) 


We get the following 

Lemma. Let t„, := max{p /jt)+3 , p /if+2 , P/„+i}; then 


T m +i < 1 




/(C/+6) 


T». 


(11) 


We can now use the quantities p f and r m to estimate the approximation 
quality: 
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Lemma. 


Pi 


(') 


fW 


< 

~ qi>) 


(12) 


Proof. Recall (2): 


x 


qW + x^ q^^ + q^"} 


Hence 


P) 


W 


/W 


< 


< 


\pf±A tl 

q( f ) + x[^q( f ') + x^2 q^*"} 9^ 

fOp?>+ x ? t <'>pP+4 , ,W ) 


(*<'>) 2 

? ( 'W' ) +*SVW ) +4VV ) 


< 


(*<'>) 2 

+ 4 W - * (, v’i 




2 pt 
~ q( f ) 


< 


Theorem. For almost all (xi, x 2 ) G B there is a constant d > 0 such that 

1 

Proof: We will first consider the case that {x\,xz) G M. By (3), (4) we 
know that 

q fi«+ 6 ) < 1 6jy ( /w+5 ) ]\[ ^ , " +4 ) ]V ( /,//+3 ) AT a/ ^ AT ('*) 

Hence 

q{fni) 1 

1 —±-< 1_ 

^(^+ 6 ) “ 16N ( / -+ 5 ) N ( / -+ 4 ) N ( / »+ 3 )AT ( / *+ 2 ) N ( / * +1 )N ('*) 


X; - ■ 


a- 


w 


rW 
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Define 


f(x u x 2 ) := log! 1 


1 


16 AT ( 6 ) N ^ NW N W N W N 0) 


Let G M; define the return time r^ii* 1 5 ^ 2 ) := 

min{>£ >0 ^(xi,x 2 ) G M}, and the induced transformation 


S M M, S m ■= S r ' ,{x " X2) 

We then have 

S'm( x i>* 2) = x 2 ). 

Since Brun’s Algorithm is ergodic and conservative (for a proof see e.g. 
Schweiger 1998), so is the system (M,Sp); we can apply the ergodic 
theorem and get 


lim - 

m— >00 m -\- 1 


m 


E 




Jiw/( y i» x 2VA* log*i n 

m(m) >(m) <u 


Thus by (11), for large enough, 

A < cKf >W) 

Since 


w-^00 m /x(M) ’ 

r„ < cK? 

and, for t large enough, 

Pt < cK\. 

We have p(M) > 0, hence (since S is conservative) /q is finite a.e., and we 
can generalize the result to (xi, x 2 ) G B. 


On the other hand, by (3) and (4) we estimate 

q {l+X) < 2N (t) q {t) < -L q {, \ 
*1 


and define 



238 B. R. Schratzberger: Exponent of Convergence for Brun’s Algorithm 


We have 


and 


Therefore 


1 


lim = : -log K 2 < o, 

'~°° * i?\ 


l {t) < K' 


Pt < 


(? (/) ) 


d 1 > 


and by Lemma (12) 


P) 


(/) 


/(') 


< 


(? w ) 


1 +d 


a.e. 
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